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ABSTRACT 

We observed PSR J0437-4715 with the FUV-MAMA detector of the Hubble 
Space Telescope Imaging Spectrometer (STIS) to measure the pulsar's spectrum 
and pulsations. For the first time, UV emission from a millisecond pulsar is 
detected. The measured flux, (2.0 ± 0.2) x 10~ 15 erg s" 1 cm" 2 in the 1150-1700 
A range, corresponds to the luminosity Lfuv = (4.7 ± 0.5) x 10 27 erg s -1 , for 
the distance of 140 pc and negligible interstellar extinction. The shape of the 
observed spectrum suggests thermal emission from the neutron star surface with 
a surprisingly high temperature of about 1 x 10 5 K, above the upper limit on 
the surface temperature of the younger "ordinary" pulsar J0108— 1431. For the 
few-Gyr-old J0437— 4715, such a temperature requires a heating mechanism to 
operate. The spectrum of J0437— 4715 shows marginal evidence of an emission 
line at 1372 A, which might be a gravitationally redshifted Zeeman component 
of the Hydrogen Lya line in a magnetic field ~ 7 x 10 8 G. No pulsations are 
detected, with a 3er upper limit of 50% on pulsed fraction. 

Subject headings: pulsars: individual (PSR J0437-4715) — stars: neutron - 
UV: stars 



1. Introduction 

Millisecond (recycled) pulsars are very old neutron stars (NSs) spun up by accretion in 
binary systems. So far, X-ray observations have been the only source of information about 
emission from millisecond pulsars (MSPs) outside the radio band (Becker & Pavlov 2001; 
Becker & Aschenbach 2002). Based on their X-ray pulse profiles and spectra, the X-ray 
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emitting MSPs can be divided into two distinct groups. The pulsars from the first group 
(e.g., PSR B1821-24, B1937+21, J0218+4232) show X-ray pulse profiles with narrow peaks, 
resembling those seen in radio, and large pulsed fractions, > 50%. They have hard power-law 
spectra, with photon indices V = 1-2, and high spin-down luminosities, E ~ 10 35 -10 36 erg 
s -1 . The X-ray radiation from these MSPs is interpreted as nonthermal emission produced 
in the pulsar magnetosphere. 

The second group consists of MSPs with smoother X-ray pulsations, lower pulsed frac- 
tions, and smaller E (~ 10 33 -10 34 erg s _1 ). In those few cases when the X-ray spectra 
are available (e.g., PSR J0437— 4715 and J0030+0451), they cannot be fitted with a single 
power-law model. The fits with simple spectral models (power-law, blackbody) require at 
least two spectral components, one of which is very soft. Likely, this soft component can be 
interpreted as thermal emission from NS polar caps, with a temperature ~ 1 MK. Such polar 
caps, heated by a backward flow of relativistic particles accelerated in the magnetosphere 
above the NS magnetic poles, are predicted by virtually all pair-cascade pulsar models (e.g., 
Ruderman & Sutherland 1975; Arons 1981; Harding & Muslimov 2002). The thermal com- 
ponent cannot be seen in the first group of MSPs because it is buried under the stronger 
nonthermal component, similar to young ordinary pulsars (e.g., Pavlov, Zavlin, & Sanwal 
2002). 

Detailed X-ray studies of PSR J0437— 4715, the brightest MSP of the second group, 
suggest that the thermal component of its radiation is emitted from a region with a nonuni- 
form temperature, decreasing from the magnetic poles towards the equator (Zavlin & Pavlov 
1998; Zavlin et al. 2002). Such a nonuniformity could be interpreted as due to a heat flow 
away from the polar cap. In addition to the external (polar-cap) heating, a variety of internal 
heating mechanisms can operate in the NS interiors, such as dissipation of the NS rotational 
energy and magnetic field (Schaab et al. 1999, and references therein). Consequently, the 
temperature distribution and, particularly, the lower value of the surface temperature depend 
on the thermal conductivity of the NS matter and the relative contributions from the exter- 
nal and internal heating. X-ray observations mainly probe thermal emission from the hot 
polar regions, being less sensitive to the emission from the rest of the NS surface with a lower 
temperature. This low-temperature emission can only be observed in the optical-UV range. 
Measuring the temperature of the NS surface in MSPs is important because it can constrain 
the NS heating models and provide information about the physical processes operating in 
the NS interior. If the magnetospheric component dominates in the optical-UV, its detection 
would help elucidate the properties of relativistic particles in the MSP magnetospheres. 

Most MSPs reside in binary systems, usually with a low-mass white dwarf (WD) com- 
panion that, as a rule, is expected to be brighter in the optical than the MSP itself. Therefore, 
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solitary MSPs look more suitable for studying the NS optical emission. However, no firm 
detections of optical/UV emission from solitary MSPs have been reported. Even very deep 
VLT observations of the solitary MSPs J0030+0451 (Koptsevich et al. 2003) and J2124-3358 
(Mignani & Becker 2003) gave negative results, putting some constraints on the nonthermal 
emission in the optical. On the other hand, if the companion of a binary MSP is sufficiently 
cold, it is very faint in the UV range, which can be used to observe NSs in nearby binary 
MSPs, particularly their thermal emission. 

The best target for such observations is PSR J0437-4715, the nearest and the brightest 
binary MSP (P = 5.76 ms, d = 139 ± 3 pc, r = P/(2P) = 4.9 Gyr, E = 4n 2 IPP- 3 = 
3.8 x 10 33 /45 er g s_1 — van Straten et al. 2001). Its binary companion is a cold WD with 
the effective surface temperature of about 4000 K (Danziger et al. 1993) and orbital period 
of 5.5 days. Optical emission from the binary is dominated by the WD (R — 20.1, V — 20.8, 
B = 22.2), making optical detection of the pulsar impossible. This prompted us to carry 
out observations of the system in the far-ultraviolet (FUV) range with the Hubble Space 
Telescope (HST). In this paper we report first detection of UV emission from a non-accreting 
MSP. The details of the observations and the data analysis are presented in §2 and §3. The 
results and their implications are discussed in §4 and summarized in §5. 



2. Observations 

PSR J0437-4715 ( J0437 hereafter) was observed on 2001 August 24 (start date 52 145.23638667 
MJD) with the Space Telescope Imaging Spectrograph (STIS). The source was imaged on 
the Far-Ultraviolet Multi Anode Micro-channel Array (FUV-MAMA). The low-resolution 
grating G140L (which covers the wavelength interval ~ 1150-1700 A) with the 52" x 0'.'5 slit 
were used. The data were taken during five consecutive orbits, including target acquisition. 
We used the WD companion of the pulsar as the acquisition target, centering it in the slit 
to a ±0'/01 nominal accuracy. The total scientific exposure time was 14,150 s (2,150 s for 
the first exposure and 3,000 s for each of the subsequent exposures). FUV-MAMA was op- 
erated in TIME- TAG mode which allows the photon arrival times to be recorded with 125 
/is resolution. 
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3. Data analysis 

For each exposure, we processed the raw, "high-resolution" images (2048 x 2048 pixels; 
plate scale of 0'/0122 per pixel — see §11 of the STIS Instrument Handbook 4 , IHB hereafter) 
using the calibration files available on January 23, 2003. As an output, we obtained flat- 
fielded "low-resolution" (1024 x 1024 pixels; plate scale / .'0244 pixel -1 ; spectral resolution 
0.58 A pixel" 1 ) images and used them for the spectral analysis. 

The processed images show a nonuniform detector background that consists of a flat 
(constant) component and the so-called "thermal glow" component (Landsman 1998) that 
dominates over most of the detector area. The thermal glow is the strongest in the upper- 
left quadrant of the detector (Fig. 1), where the dark count rate can exceed the nominal 
value, 6 x 10~ 6 counts s -1 pixel -1 , by a factor of 20. The brightness of the thermal glow 
increases exponentially with increasing the FUV-MAMA power supply temperature above 
38.9 K (Landsman 1998). In our observation, this temperature was rising from 40.89 K 
during the first orbit to 42.36 K during the last orbit, resulting in brightening thermal glow. 
In addition, the overall shape and small-scale structure of the thermal glow varied slightly 
between the exposures. 

Because of the high background, the dispersed source spectrum is hardly discernible in 
the separate exposures. Nevertheless, we find the spectrum centered at Y = 351 ± 2 pixels 
in each of the flat-fielded images (the centroid slightly varies with X), where X and Y arc 
the image coordinates along the dispersion and spatial axes, respectively. This is within 
2 pixels (C/OS) of the nominal position where the target was expected to be acquired. At 
this location on the detector, the contribution of the thermal glow to the background is still 
quite large — the background averaged along the dispersion direction exceeds the nominal 
background by a factor of 10. To improve the signal-to-noise ratio (S/N), we combined the 
images from five exposures into a single image using the STSDAS 5 task "mscombine". The 
V-positions of the centroids differ by less than 3 pixels for different exposures and different 
wavelengths (X-positions). The spectrum of the source is clearly seen in the summed image 
shown in Figure 1. 

Accurate subtraction of the enhanced, nonuniform background is crucial to measure the 
spectrum of our faint target. The spectral extraction algorithm implemented in the standard 
STIS pipeline (task X1D) does not adequately correct for the nonuniform background while 
extracting the spectrum of such a faint source. Therefore, we developed an IDL routine with 
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additional capabilities of grouping and fitting the background to reduce the data. 

Since the source spectrum occupies only a small region on the detector, we do not 
attempt to subtract the background globally. Instead, we scan the count distribution within 
two strips, 324 <Y< 343 and 362 < Y < 381, adjacent to the source region, 344 < Y < 361. 
To obtain the spectrum with a sufficiently high S/N, we have to bin the spectrum heavily; 
after some experimenting, we chose four spectral bins (A-bins): 1155-1187, 1248-1283, 1316— 
1376, and 1400-1702 A. The first two A-bins are chosen to avoid contamination by the 
geocoronal line emission (Lya line at 1216 A and the Oxygen line at 1304 A; the other 
geocoronal Oxygen line, at 1356 A, is not seen in this observation). Because of an enhanced 
background at 1376-1400 A, we also exclude this region from the spectral analysis, which 
determines the choice of the third A-bin (1316-1376 A). The remaining counts are grouped 
into a single bin (1400-1702 A) to have comparable S/N in the second through fourth A-bins. 

For each of the A-bins, we calculate the total number of counts, N t , within the extraction 
boxes of different heights (one-dimensional apertures): A s = 3, 5, 7, 9, 11, 13, 15, and 17 
pixels, centered at Y — 351, 352, 353, and 352 for the 1st, 2nd, 3rd, and 4th A-bins, 
respectively. To evaluate the background, we first clean the background strips (see above) 
from outstanding (> 10~ 3 cts s -1 pixel -1 ) values ("bad pixels") by setting them to local 
average values (median filtering or standard sigma-clipping algorithms are not applicable in 
this case — see IHB §7.4.2). Then, for each of the A-bins, we fit the ^-distribution of the 
background counts with a first-order polynomial (interpolating across the source region), 
estimate the number iV b of background counts within the source extraction aperture A s , and 
evaluate the number of source counts, N s = N t — Nb (Table 1). 

The uncertainty 5N S of the source counts can be evaluated as 5N S = [N t + (SNb) 2 ] 1 ^ 2 , 
where SNb is the background uncertainty in the source aperture. For a uniform, Poissonian 
background, this quantity can be estimated as = [Nb(A s /A b )] 1/2 , where A b (= 40 

pixels) is the aperture where the background was measured. To account for the background 
nouniformity, we estimated 5Nb directly from the image. We binned the distribution of 
background counts along the F-axis with the bin sizes equal to A s (for A s = 3, 5, and 
7) and calculated 5N b as root-mean-square of the differences between the actual numbers 
of background counts in the bins and those obtained from the fit to the background. The 
background uncertainties obtained in this way are systematically larger (see Table 1) than 
the uncertainties estimated assuming a uniform, Poissonian background. 

We find that the signal-to-noise ratio, S/N = N S /5N S , calculated for each A-bin and 
for different extraction box heights, is the largest for A s = 3 or 5 pixels, depending on the 
A-bin. Since the results at small heights are sensitive to the deviation of the spectrum from 
the straight line (along the dispersion direction) and possible small misalignments between 
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the frames taken in the different orbits, we choose A s = 5 pixels as the optimal value. We 
adopt this source extraction aperture for further analysis, which contains 53%, 61%, 64%, 
and 65% of source counts for the first through the fourth A-bins, respectively). 

For each A-bin, we calculate the average spectral flux defined as 

f AX R x \F x d\ _ Cl 
{ Xh ~ U^AdA = / AAi i?AAdA' (1) 

where Cj is the source count rate in the 2-th A-bin, and R\ is the system response (which 
includes the Optical Telescope Assembly throughput and accounts for the grating and slit 
losses and the finite size of the source extraction aperture; see §3.4.12 of the HST Data 
Handbook for STIS 6 for details). The resulting flux values are given in Table 1, while the 
spectrum is shown in Figure 2. The total flux in the 1150-1700 A range (AA = 550 A), 
can be estimated as F ~ AA (E^^AA;) (£\ AA;) -1 ~ (2.0 ± 0.2) x 10~ 15 erg s" 1 cm" 2 , 
corresponding to the luminosity L FUV = AncPF = (4.7 ± 0.5) x lQ 27 d\ i0 erg s _1 . 

In the original spectrum we see a brightening around 1372 A (marked with a vertical 
arrow in Fig. 1), resembling an emission line. The count rate spectrum of the possible line 
and its vicinity, binned to 1.16 A (4 high-resolution pixels) bins, is shown in Figure 3. The 
count excess, 30 ± 10 counts in the three bins, corresponds to the flux (4.0 ± 1.3) x 10~ 17 
erg cm~ 2 s _1 in the 1370.1-1373.6 A range. This spectral feature, detected at a 3 er level, 
is clearly not of a geocoronal origin (the nearest geocoronal oxygen line, at 1356 A, is not 
seen in this observation), but we cannot rule out the opportunity that it may be associated 
with an anomaly in the FUV-MAMA background seen right above the feature (in the cross- 
dispersion direction). Therefore, we consider this as a marginal detection until the feature 
is confirmed in another observation. 

For the timing analysis, we used the so-called TIME- TAG data files that contain the 
photon arrival times, recorded at a 125 /is time resolution, and high-resolution detector 
coordinates (see §2) associated with each of the events. We extracted 2964 events from 
the above-defined four A-bins, with the height of extraction box equal to 9 high-resolution 
pixels (~ 22.5% of these counts are expected to come from the source). Note that the 
background count rate in these data is about 3% higher than that in the data used for 
the spectral analysis where "outstanding pixels" were filtered out (see above). The arrival 
times of the events were corrected for the Earth and spacecraft motions and transformed to 
barycentric dynamical times (TDB) at the solar system barycenter, using the STSDAS task 
"odelaytime" . To correct the TDB arrival times for the effect of the binary orbital Doppler 
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shift (Taylor & Weisberg 1989) and to extract the light curve, we used the pulsar and binary- 
orbit ephemerides from van Straten et al. (2001). 

Since the frequency of radio pulsations, / = 173.687948857032 Hz at the epoch of the 
STIS observation, is known with high precision, much higher than one can achieve in the 
time span T span = 24 638 s of our observation, there is no need to search for the period in the 
FUV data. We folded the corrected arrival times with the pulsar's period and obtained the 
light curves for different numbers of phase bins. An example of the light curve for 4 phase 
bins, with an arbitrarily fixed reference phase (histogram) and averaged over the reference 
phase (smooth line; see Gregory & Loredo 1992 for a description of this procedure) is shown 
in Figure 4. The observed pulsed fraction can be estimated as f p = 4^%, which means that 
the pulsations are not statistically significant. 

The upper limit on the pulsed fraction can also be estimated from the Z\ test (e.g., 
Buccheri et al. 1983) that gives Z 2 n = 2.7, 2.8, 5.2, 8.1, and 8.7 for n = 1, 2, 3, 4, and 5, 
respectively, at the radio pulsation frequency {n is the number of harmonics included). For 
randomly distributed arrival times, the probabilities of getting Z 2 n larger than these values 
are 0.26, 0.59, 0.52, 0.42, and 0.56, respectively, which confirms that no periodic signal is 
detected. An upper limit on the pulsed fraction can be estimated from the Z\ value, assuming 
sine-like pulsations. At a 3<r confidence level, the limit is 11%— 12%, depending on the method 
used for estimation (e.g., Protheroe 1987; Vaughan et al. 1994; Brazier 1994). It translates 
into an upper limit /p nt < 48%-53% for the intrinsic pulsed fraction of J0437, corrected for 
the background contribution. The upper limit is close to the pulsed fraction at E > 1 keV, 
and it exceeds the pulsed fraction, ~ 30%, observed in soft X-rays (Zavlin & Pavlov 1998). 
If the FUV radiation is thermal, as suggested by the spectrum (see §4), we expect the pulsed 
fraction to be even lower than 30% (because the FUV radiation is presumably emitted from 
a large fraction of the NS surface), so the limit is not truly restrictive. The reason for that is 
the very high detector background. We attempted to improve the S /N repeating the analysis 
for the first four orbits, with the lower detector background, but obtained approximately the 
same upper limit. 



4. Discussion 

4.1. Neutron star or white dwarf? 

First of all, we should understand whether the observed FUV spectrum is emitted by 
the WD companion or by the NS. The results of the optical photometry (Danziger et al. 
1993; Bailyn 1993) yield T cS « 3750-4500 K for the WD. The dashed line in Figure 2 
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shows a blackbody spectrum (for T = 4000 K, R = O.O25_R , d = 140 pc) that crudely 
fits the photometric data and clearly falls well below the FUV data points. More accurate 
comparison should use cool WD atmosphere spectra (e.g., Bergeron et al. 1995; Hansen 
1999), which can be harder than the blackbody at UV frequencies. We have applied a 
number of such models (provided by P. Bergeron 2003, private communication) for various 
chemical compositions, effective temperatures and gravitational accelerations, scaling them 
to fit the optical data, and found that the observed FUV flux exceeds the WD atmosphere 
model predictions by several orders of magnitude in all the cases. Thus, we conclude that 
the observed FUV emission originates from the NS. 



4.2. Thermal origin of the FUV spectrum 

We fit the spectrum with the absorbed power-law model, Fa = Fi 50 o (A/1500 A)" A x 
Iq~oaa(\)e(b-v)^ w h ere the ultraviolet extinction curve A(\) is adopted from Seaton (1979). 
The (low) interstellar extinction in the direction of J0437 is poorly known; different authors 
use the color excess E(B — V) values from < 0.01 (Bell et al. 1993) to 0.07 (Danziger et 
al. 1993). We performed the fits for three values, E(B - V) = 0, 0.03, and 0.07, and found 
the power-law indices ot\ = —4.0 ± 1.2, —4.2 ± 1.2, and —4.4 ± 1.2, and the normalizations 
F 1500 = 2.6 ±0.4, 3.4 ±0.5, and 4.4 ± 0.7 x 10~ 18 erg cm" 2 s" 1 A" 1 , respectively (see Fig. 5). 
The inferred slope ot\ (albeit rather uncertain because of the strong detector background) 
is close to that of the Rayleigh- Jeans spectrum, F\ oc A~ 4 . It indicates that, most likely, 
the observed radiation is thermal radiation from the NS surface rather than magnetospheric 
radiation. At d = 140 pc, the normalization parameter of the Rayleigh- Jeans fit corresponds 
to the brightness temperature T 5 = 0.215 F^oi?^ 2 K , where T 5 = T/(10 5 K), R 13 is the 
radius of the emitting sphere in units of 13 km. (Here and below, the temperature and 
the radius are given as measured by a distant observer.) Our fits (with fixed ct\ = —4) give 
T 5 Rj 3 = 0.62±0.06, 0.79±0.08, and 1.1 ±0.1, for E(B-V) = 0, 0.03, and 0.07, respectively. 

At T < 10 5 K the blackbody spectrum deviates appreciably from the Rayleigh- Jeans 
limit in the FUV range. To investigate the range of lower temperatures, we also fit the 
absorbed blackbody model to the observed spectrum. The confidence contours in the T- 
R plane (Fig. 6) show that the lower limit on the surface temperature is 4.2 x 10 4 K, 
corresponding to a radius of 37 km, at a 68% confidence level (3.1 x 10 4 K and 59 km at a 90% 
level). For a typical NS radius R = 13 km, the inferred surface temperature is 0.107± 0.010, 
0.121 ± 0.012, and 0.156 ± 0.014 MK, for E(B - V) = 0, 0.03, and 0.07, respectively. The 
corresponding bolometric luminosity can be estimated as L bo i = 1-2 x 10 29 T 5 4 i? 2 3 erg s _1 ; 
for instance, L hoi = (2.6 ± 1.0) x 10 29 erg s" 1 for E(B -V) = 0.03, R = 13 km. 
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Because of the large uncertainties of the FUV data points, deviations of the source 
spectrum from the Rayleigh- Jeans limit cannot be established from the FUV data alone. 
Therefore, to infer the upper limit on T (lower limit on R), we have to use the X-ray data. 
For the single-temperature blackbody model, we find that T should be lower than 0.23 
MK (R > 8.0 km), for E(B — V) = 0.03, in order not to exceed the observed soft X-ray 
spectral flux (see Fig. 7). It is clear from Figure 7 that the FUV and X-ray data cannot 
be described by a single-temperature model. Based on the X-ray observations, Zavlin et al. 
(2002) suggested a two-temperature model: T core = 2.1 MK, R com = 0.12 km; T rim = 0.54 
MK, i? rim = 2.0 km, where the "core" and the "rim" correspond to central and peripheral 
parts of the pulsar's polar caps. We see from Figure 7 that extension of this model to the 
FUV range underestimates the observed FUV spectral flux. Most probably, this means that 
the NS surface temperature is a function of the polar angle, such that it decreases from ~ 2 
MK at the magnetic poles to ~ 0.1 MK at the bulk of the surface. Of course, the estimates 
for the brightness temperature above do not imply a uniform surface temperature away 
from the polar caps. To determine the dependence of the surface temperature on magnetic 
colatitude, phase-resolved spectroscopy is required. 



4.3. Internal heating 

The inferred surface temperature, T ~ 10 5 K, is much higher than the temperatures 
< 10 3 K expected for a passively cooling few-Gyr-old neutron star (e.g., Tsuruta 1998; 
Schaab et al. 1999 [S99 hereafter]). Therefore, it requires a heating mechanism to operate 
in J0437. If the sources of heat are in the highly conductive (hence almost isothermal) 
NS interiors, the surface temperature of 0.1 MK implies the interior temperature of order 
1 MK (Gudmundsson, Pethick, & Epstein 1983). Various mechanisms of internal NS heating 
were investigated in a number of works (e.g., S99, and reference therein). The mechanisms 
relevant to MSP heating can be divided into two groups. Firstly, heat can be produced by 
the dissipation of energy of differential rotation caused by frictional interaction between the 
faster rotating superfluid core and the slower rotating outer solid crust (Shibazaki & Lamb 
1989; Larson & Link 1999, and references therein). The heating mechanisms from the second 
group are associated with readjustment of the NS structure to a new equilibrium state as 
star rotation slows down. If the NS crust is solid, it will undergo cracking when the tension 
force exceeds a critical value, and the crust strain energy will be released as heat (Cheng 
et al. 1992). In addition, as the star spins down, the centrifugal force decreases and central 
density increases causing the shift in chemical equilibrium, which modifies the rate of nuclear 
reactions and may lead to heat release (Reisenegger 1995). 
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The high surface temperature of J0437 strongly constrains heating mechanisms. For 
instance, such a temperature rules out crust cracking and chemical heating of the crust as 
the main heat sources (cf. Fig. 7 of S99). Chemical heating of the core and frictional heating 
remain viable mechanisms. For the latter, one can constrain the excess angular momentum 
AJ S , residing in the superfluid, and the angular velocity lag lD, between the superfluid and 
the crust, averaged over superfluid moment of inertia. Using equation (9) of Larson & Link 
(1999), we obtain AJ S = L hol /\tt\ = 1.1 x 10 41 T£R 2 13 erg s, cu = AJ S /I S = 0.157£R? 3 rad 
s _1 , where (l = —1.086 x 10~ 14 rad s -2 is the time derivative of the angular frequency of the 
pulsar, I s = 7.3 x 10 43 g cm 2 is the moment of inertia of the differentially rotating portion 
of the superfluid, estimated for the Friedman & Pandharipande (1981) equation of state. 
Because other mechanisms can contribute to heating of J0437, these estimates of AJ S and 
uj should be considered as upper limits. 



4.4. External heating 

In addition to the above-discussed internal sources, heat can be provided by relativis- 
ts particles created in the pulsar's acceleration zones and bombarding the NS surface 
(e.g., Harding & Muslimov 2002, and references therein). The energy of these particles 
is released in polar caps at the NS magnetic poles, with a radius r pc ~ (R 3 tt/ 'c) 1 / 2 ~ 
2(i?/10 km) 3 / 2 (P/6 ms)" 1 / 2 km. Moreover, inward-directed radiation from this particle 
population or even precipitating secondary particles created in the closed field line zone 
can heat the surface outside of the traditional polar caps (Wang et al. 1998). The energy 
dissipated in the NS polar regions per unit time (heating luminosity) depends on the ge- 
ometry of the acceleration region and the degree of electric field screening. The numerical 
calculations, carried out by Harding & Muslimov (2002) for the case of MSPs, give heating 
rates somewhat lower than (but close to) the bolometric luminosity inferred from the X-ray 
observations of J0437, L^ c ol ps 2 x 10 30 erg s~ 4 (Zavlin et al. 2002). A major fraction of the 
thermal energy deposited by the relativistic particles propagates towards the surface and is 
radiated in soft X-rays from the polar cap(s). However, some energy flows inwards, heating 
deeper NS layers beneath a larger surface area. If a fraction of this energy reaches the inner 
crust, it spreads over the whole interior of the NS because of its very high thermal conduc- 
tivity, so eventually it will be radiated from the entire NS surface, together with the heat 
supplied by possible internal heating mechanisms. The effect of the external heating on the 
temperature distribution in an NS requires numerical solution of a complicated problem of 
heat transport, with allowance for radiation from the surface. This problem has not been 
solved, to the best of our knowledge, so it is hard to assess the contribution of the polar 
cap heating to the observed FUV radiation. We can only state that if the external heating 
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dominates, then a fraction ~ / L^ c ol ~ 0.1-0.5 of the energy deposited at the polar cap 
goes to heating of the NS. Crude estimates show, however, that such a large fraction of heat 
flowing inwards would require thermal conductivities much higher than calculated for the 
NS crust (e.g., Jones 1978), suggesting that particle precipitation only on the open field line 
zone cannot account for general surface heating of old NSs. 

In principle, a very old NS could be heated by exotic particles (or products of their 
decay) captured in the NS matter. For instance, Hannestad, Keranen, & Sannino (2002) 
consider NS heating by products of decay of Kaluza-Klein majorons and gravitions, that 
form a halo around the NS, and conclude that the thermal luminosity of an old NS should 
reach a constant level corresponding to the energy deposited by decay particles. (See Pavlov, 
Stringfellow, & Cordova 1996 for references to earlier papers.) Although the Occam razor 
suggests that less exotic mechanisms should be explored first, the possibility to indirectly 
investigate exotic particles through thermal radiation from MSPs cannot be dismissed at the 
present stage. 

4.5. Are millisecond pulsars hotter than ordinary old neutron stars? 

It is interesting to compare the inferred surface temperature of J0437 with the tempera- 
tures (or upper limits) for other old pulsars, including MSPs. We compiled the results of the 
available UV-optical observations of old (> 1 Myr), nearby (< 300 pc) NSs in Table 2 and 
re-estimated the temperatures using the most recent results on the distances. In addition to 
the old pulsars, we included the famous "dim" isolated NS, RX J1856.5— 3754, whose age is 
unknown, because it might be an MSP (Pavlov & Zavlin 2003). We see from Table 2 that 
the temperature of J0437 is the lowest measured temperature for a NS. The observations 
of the other two nearby MSPs, J2124— 3358 and J0030+0451, were not nearly as deep as 
required to detect their thermal radiation at a similar temperature. All the upper limits 
on T in Table 2 are also above the temperature of J0437 with one exception — T 5 < 0.88 
for PSR J0108— 1431, whose spin-down age is a factor of 30 smaller than that of J0437 3 . 
Although the difference in the temperatures is not very large, it should be noted that the 
upper limit in Table 2 (see also Mignani et al. 2003b) is rather conservative, in the sense that 
it is obtained assuming d = 200 pc, the largest among the distances, 60-200 pc, estimated 
from various models of Galactic electron distribution (Tauris et al. 1994; Taylor & Cordes 



3 Although the true age of a pulsar may differ substantially from its spin-down age, the low temperature 
of the WD companion implies (e.g., Hansen & Phinney 1998) that J0437 is older than 3 Gyr, that is a factor 
of 18 larger than the spin-down age of PSR J0108-1431. 
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1993; Cordes & Lazio 2002). For instance, the limit is T 5 < 0.45 for d = 130 pc (Mignani et 
al. 2003a). Even if the true distance is somewhat larger than 200 pc, the upper limit remains 
surprisingly low, given the much younger spin-down age of JO 108— 1431. 

The higher temperature of the older J0437 cannot be attributed to accretion in the 
binary system because the accretion ceased about 3 Gyrs ago, while the time for a NS to 
cool to 10 5 K, losing the accretion-generated heat, is < 10 Myr. The high temperature can 
possibly be explained by stronger frictional heating, consistent with a factor of 15 larger |f2|. 
It is also possible that "rotochemical heating" of the NS core (Reisenegger 1995,1997) has 
been turned on in this MSP, or J0437 is simply sufficiently old to accumulate enough heat 
by capturing exotic particles (see §4.4). Finally, the higher temperature of J0437 might be 
associated with stronger magnetospheric heating, since its E is about 700 times larger than 
that of J0108— 1431. In this case one should expect MSPs with higher E (up to 3 orders of 
magnitude in most powerful MSPs, such as PSR B1821— 24) to be even hotter than J0437, 
but it would be difficult to discern their surface emission since these pulsars are also strong 
sources of magnetospheric radiation. 

Whatever is the reason of the higher temperature of J0437, it is tempting to assume 
that the surface temperatures of MSPs are generally higher than those of old "ordinary" 
pulsars of similar or even younger ages. Since ordinary isolated NSs cool down with age even 
in the presence of internal heating (Tsuruta 1998; S99), this assumption implies that the 
transition from ordinary pulsars to MSPs is accompanied by a considerable growth of NS 
temperature. Moreover, the temperature is maintained high in the course of the Gyr-long 
MSP thermal evolution by a heating process, probably associated with specific properties 
of MSPs, such as fast rotation or old age. To understand the apparently "nonmonotonous" 
thermal evolution of old NSs, thermal emission from a larger sample of these objects should 
be investigated. 

4.6. Possible spectral line at 1372 A 

Although the detection of the emission line at 1372 A (hu = 9.04 eV) is marginal, it is 
interesting to speculate as to about its origin. Firstly, this could be an electron cyclotron 
line (e.g., formed in a corona above the polar cap) in the magnetic field B = 7.8 x 10 8 (1 + z) 
G, where z is the gravitational redshift. Such a field is close to the "conventional" magnetic 
field, B p = 6.6 x 10 8 Rq 3 ll[ 2 (sin a)^ 1 G, at the magnetic poles of an NS losing energy via 
magnetodipole radiation (R = 10 6 i?6 cm and / = 10 45 /45 g cm 2 are the NS radius and 
moment of inertia, a is the angle between the magnetic and rotational axes). The main 
difficulty with this interpretation is the small width of the observed line, ~ 3 A. First, 
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it implies a very uniform field in (hence, a small size of) the emitting region, AB/B < 
AA/A. For instance, if the emitting region is a hot spot at the magnetic pole, its size should 
be smaller than ~ 0.5 km. Second, the thermal (Doppler) width of the cyclotron line is 
AAb = A(2/cT/m e c 2 ) 1 / 2 | cos cos#| A, where 9 is the angle between the line of 

sight and the direction of the magnetic field. For the thermal width to be smaller than 
AA, a phase-averaged | cos#| must be < 0.12 T 6 [9 > 83° for T 6 = 1), which is hard to 
reconcile with the constraints on the rotation and magnetic inclinations obtained from the 
radio polarimetry (e.g., Manchester & Johnston 1995). 

Alternatively, if the NS surface is covered with a hydrogen atmosphere with overheated 
outer layers (e.g., due to convection), the observed line may be one of the three gravitationally 
redshifted Zeeman components of the Lya line in a magnetic field B ~ 10 8 -10 9 G. For 
instance, at B = 7 x 10 8 G, the wavelengths of the Zeeman components at the NS surface 
are 730, 1069, and 1334 A (Ruder et al. 1994). If the observed line is the redshifted vr- 
component, the redshift is z = 0.28, so the wavelengths of the redshifted a-components, 
937 and 1712 A, are outside of the observed FUV range. The thermal width of the Zeeman 

1/2 > 

component, ~ 0.6 T 6 A, is much smaller than the observed width. The line broadening 
can be caused by a magnetic field nonuniformity: AA = 3 A corresponds to AB/B ~ 0.3 
(estimated from Fig. A. 2.1 of Ruder et al. 1994), that is the putative corona may spread 
over a substantial fraction of the NS surface. The luminosity in the line, L 1372 ~ 0.9 x 10 26 
erg s -1 , is a small fraction, 2%, of the observed FUV luminosity. If the emitting region is 
in collisional (coronal) equilibrium, the corresponding emission measure is n 2 H V ~ 10 49 -10 50 
cm -3 for T ~ 0.1-1 MK. If the 1372 A line is confirmed in a future FUV observation, the 
most convincing confirmation of the Lja interpretaion would be detection of another Zeeman 
component (e.g., the a-component in the NUV range), which would allow one to measure 
both the magnetic field and the gravitational redshift. 



4.7. Magnetospheric component 

Zavlin et al. (2002) found that, in addition to the two-temperature thermal component, 
a power-law component with a photon index T = 2.2^Qg is needed to fit the X-ray spectrum 
at energies above 2 keV (Fig. 7). However, extension of this component with the best-fit 
T to the FUV range is well above the observed FUV spectrum, and even for T = 1.6 the 
extension is only marginally consistent with the FUV data (Fig. 7). We can crudely estimate 
the upper limit on the magnetospheric luminosity is the FUV range as -^puv ~ 4.2 x 10 27 
erg s _1 . The apparent discrepancy can be explained assuming that the spectrum of the 
magnetospheric radiation breaks down when the frequency decreases from X-rays to the 
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FUV range. The optical upper limits for MSPs J0030+0451 (Koptsevich et al. 2003) and 
J2124— 3358 (Mignani & Becker 2003), well below the extensions of their X-ray power- law 
components, suggest that it might be a common property of MSPs, contrary to, e.g, ordinary 
middle-aged pulsars. Alternatively, as Zavlin et al. (2002) noticed, it is possible that the 
excess counts at higher energies in MSP spectra might be interpreted as thermal radiation 
from a very hot, small area within the polar cap (e.g., T ~ 12-15 MK, R ~ 1-2 m for J0437). 
To understand which of the two interpretations is correct, deeper observations of J0437 in 
hard X-rays are required. 



5. Summary and conclusions 

The STIS/FUV-MAMA observation of J0437 provided first firm detection of an MSP in 
the optical-UV range. The FUV spectrum is best interpreted as thermal emission from the 
NS surface with a temperature of about 0.1 MK. This temperature exceeds the upper limit on 
the temperature of the younger, but less luminous, ordinary pulsar JO 108-1431. This is likely 
associated with a difference in spindown-driven heating. If magnetospheric heating plays a 
role, it must be effectively communicated, perhaps by radiation or secondary particles, to 
the bulk of the NS surface. Evolutionary differences between ordinary pulsars and MSPs 
might plausibly affect the internal thermal history. To understand thermal evolution of old 
NSs, more MSPs and ordinary old pulsars should be observed in the optical-UV range. 

Comparison of the FUV and X-ray spectra shows that the temperature is not uniformly 
distributed over the NS surface. The X-ray observations, sensitive to higher temperatures, 
show a smaller size of the emitting region, naturally interpreted as a pulsar polar cap, 
perhaps also with a nonuniform temperature. To understand the temperature distribution 
over the NS surface, phase-resolved spectroscopy in both X-rays and FUV is needed. We 
failed to detect FUV pulsations because the source was placed at a region of high detector 
background. An FUV observation of J0437 with an optimal positioning on the detector could 
detect pulsations (or put a stringent limit on the pulsed fraction) and provide information 
on the temperature distribution. 

The FUV upper limit on the nonthermal (magnetospheric) component, observed in hard 
X-rays, suggests a spectral turnover of this component at EUV wavelengths, which can be 
a generic property of MSPs. However, the upper limit is not very strong because of large 
errors associated with the high detector background. To tightly constrain the nonthermal 
component, another FUV-MAMA observation with improved S/N as well as deep NUV and 
X-ray exposures are required. 
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The marginally detected emission line at 1372 A can be interpreted as an electron 
cyclotron line or, more likely, a Zeeman component of the Hydrogen Lya line in a magnetic 
field of ~ 10 9 G. Confirming this line would be of profound importance as it provides an 
opportunity to directly measure the MSP magnetic field and gravitational redshift (if another 
Zeeman component is also detected). 
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Table 1: Counts and fluxes in A-bins 



A-bin 


N t 


iV 6 a 


SN b 


5N£ ni 


N s 


SN S 


S/N 


F\ h 


(A) 


















1155-1187 


367 


313 


9.3 


6.2 


54.3 


21.3 


2.5 


7.5 ±2.6 


1248-1283 


615 


465 


15.0 


7.6 


149.3 


28.0 


5.1 


5.4 ±0.9 


1316-1376 


816 


609 


22.6 


8.7 


206.7 


36.5 


5.8 


4.6 ±0.6 


1400-1702 


1369 


1078 


25.1 


11.6 


290.6 


44.7 


6.5 


2.9 ±0.4 


Summed 


3167 


2466 


38.2 


17.5 


700.8 


68.0 


10.3 


3.7 ±0.3 



a Number of background counts scaled to the 5 pixel extraction box. 

h Observed average spectral flux, in units of 10~ 18 erg s -1 cm~ 2 A" 1 , corrected for the finite aperture. 
c Values for summed A-bins. 
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Table 2: Blackbody temperatures of nearby old neutron stars 



NS 




Distance 15 


Age c 


F d 


A c 


E(B — 


V) { Ref. g 




(10 5 K) 


(PC) 


(Myr) 


(10" 31 cgs) 


(/iin) 






B1929+10 


< 4.9 


331 ± 10 


3.2 


5.4 


0.24 


0.10 


1 


B0950+08 


< 2.5 


262 ±5 


16 


5.1 


0.31 


0.01 


2 


J0108-1431 


< 0.88 


200 


170 


< 1.5 


0.49 


0.01 


3 


J1856-3754 


5.6 h 


117 ± 12 


? 


240 


0.15 


0.03 


4 


J0437-4715 


1.2 


139 ±3 


4900 


20 


0.15 


0.03 


5 


J2124-3358 


< 4.6 


270 


7200 


< 3.3 


0.49 


0.05 


6 


J0030+0451 


< 9.2 


320 


7700 


< 4.8 


0.49 


0.05 


7 



"Blackbody (brightness) temperature from optical/UV observations, for an NS radius of 13 km and quoted 
distance and color index. 

b Distances estimated from parallaxes (with ± uncertainties; see Brisken et al. 2002; Walter & Lattimer 2002; 
van Straten et al. 2001) or from the dispersion measure (Cordes & Lazio 2002). 
c Characteristic age of pulsars, t = P/{2P). 

^Observed spectral flux (or upper limit), in units of 10~ 2 /ttJy, at the wavelength quoted in the next column. 
e Rcfcrence wavelength. 
•^Adopted color index. 

^References for the spectral flux measurements. -(1) Mignani et al. 2002; (2) Pavlov et al. 1996; (3) Mignani 
et al. 2003b; (4) Pons et al. 2002; (5) this work; (6) Mignani & Becker 2003; (7) Koptsevich et al. 2003. 

''Lower temperatures, T 5 < 3.9, corresponding to R > 16 km, are required in a two-component blackbody 
model (Pavlov et al. 2002; Brajc & Romani 2002). 
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Fig. 1. — Distribution of counts on the FUV-MAMA detector. The spatial (Y) and dispersion 
(X) axes are in the vertical and horizontal directions, respectively. The spectrum of J0437 
is shown by horizontal arrows. The background is clearly dominated by the nonuniformly 
distributed "thermal glow" which is the strongest at the upper left corner and the weakest 
at the bottom. The vertical arrow shows the location of the possible spectral feature (see 
§4.6 and Fig. 3). 
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Fig. 2. — FUV spectrum of J0437. The error bars on the left represent the measured average 
fluxes in four A-bins. The dash-dotted line is the fit with the absorbed power-law model for 
E(B — V) = 0.03. The dotted line shows the same model but dereddened. The dashed line 
is a blackbody spectrum with the temperature of 4000 K, fitting the dereddened B,V,R,I 
WD fluxes, taken from Danzinger et al. (1993). 
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Fig. 3. — Total (source + background) count rate distribution in the region around the 
possible emission feature at 1372 A. 
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Fig. 4. — Pulsar light curve for 4 phase bins. The smooth line shows the light curve averaged 
over the reference phase (see §2). 
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Fig. 5. — Confidence levels (68% and 90%) for the absorbed power-law model fit, for E(B — 
V) = 0, 0.03, 0.07. 




Fig. 6. — Confidence levels (68% and 90%) for the absorbed blackbody model fit, for E(B — 
V) = 0.03. The contours are cut (the dashed line in the top left corner) because the X-ray 
model flux exceeds the observed one at higher temperatures and smaller radii (see Fig. 7). 
The horizontal dashed line shows the temperature corresponding to a plausible NS radius of 
13 km. 



-26- 



B 

o 



0) 



El, 
GO 

o 



-12.0 



-12.5 - 



-13.0 - 



-13.5 



-14.0 - 



i — i — i — i — i — i — i — i- 



i 

- \ 
\ 
i 
i 
i 
i 



/ 
/ 

7 / 

FUV^ " " >' 




41 / \....- 



\ 

i0.12 MK 1 



\ 



- 14.5 h jx--'fi: 

»•" I / /' 
1 1 I' 



\ 0.042 MK 



\ 



13 km 



\ 37 km \ 
\ 



-i i i i i i i i_ 



15.0 15.5 16.0 16.5 17.0 

Log Frequency, Hz 



17.5 



18.0 



Fig. 7. — Multiwavelength spectrum of J0437. The data points on the left represent the 
measured FUV fluxes; the unabsorbed (dereddened) fluxes are only slightly higher for plau- 
sible values of E(B — V) (see Fig. 2). The data points on the upper right are the unabsorbed 
X-ray fluxes from Chandra ACIS and ROSAT PSPC observations (Zavlin et al. 2002). Three 
dashed lines marked with temperature and radius values show the range of formally accept- 
able blackbody models (see §4.2 for details). The dashed line marked "WD" corresponds to 
the blackbody with the temperature of 4000 K (cf. Fig. 2). The dash-dotted line that goes 
through the X-ray data points is the two-temperature hydrogen atmosphere + power-law 
model (see §4.2, §4.7, and Zavlin et al. 2002). The power-law component with the photon 
index T = 1.6 is shown separately by the dotted line. 



